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Abstract 

, Electromagnetic wave scattering by many parallel to z— axis, thin, 

^ [ impedance, circular infinite cylinders is studied asymptotically as a — > 0. 

Let D m be the crossection of the m—th cylinder, a be its radius, and 
&m = (x m i,x m 2) be its center, 1 < m < M, M = M(a). It is assumed 
that the points x are distributed so that 



A/-(A) = ^i-^ N(x)dx[l + o(l)], 



where A/"(A) is the number of points x m in an arbitrary open subset A 
of the plane xoy. The function N(x) > is a given continuous function. 
An equation for the self- consistent (efficient) field is derived as a — > 0. 
A formula is derived for the effective refraction coefficient in the medium 
in which many thin impedance cylinders are distributed. These cylinders 
may model nanowires embedded in the medium. Our result shows how 
these cylinders influence the refraction coefficient of the medium. 

PACS: 03.40.Kf, 03.50 De, 41.20.Jb, 71.36,+c 
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1 Introduction 

There is a large literature on electromagnetic (EM) wave scattering by an array 
of parallel cylinders (see, e.g., [2J, where there are many references given, and 
[3]). Electromagnetic wave scattering by many parallel to z— axis, thin, circular, 
of radius a, infinite cylinders, on the boundary of which an impedance boundary 
condition holds, is studied in this paper asymptotically as a — > 0. The cylinders 
are thin in the sense ka <IC 1, where k is the wave number in the exterior of the 
cylinders, 

The novel points in this paper include: 

1) The solution to the wave scattering problem is considered in the limit 
a — > when the number M = M(a) of the cylinders tends to infinity at a 
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suitable rate. The equation for the limiting (as a — > 0) effective (self-consistent) 
field in the medium is derived, 

2) This theory is a basis for a method for changing refraction coefficient in 
a medium. The thin cylinders model nanowires embedded in the medium. The 
basic physical resul of this paper is formula (|42[) , which shows how the embedded 
thin cylinders change the refraction coefficient n 2 {x). 

Some extension of the author's results ([7]-[13]) is obtained for EM wave 
scattering by many thin perfectly conducting cylinders. The techniques used 
are similar to the ones developed in |14) . 

Let D m , 1 < m < M, be a set of non- intersecting domains on a plane P, 
which is xoy plane. Let x m £ D m , x m — (x m i, x m 2), be a point inside D m and 
C m be the cylinder with the cross-section D m and the axis, parallel to z-axis, 
passing through x m . We assume that x m is the center of the disc D m if D m is 
a disc of radius a. 

Let us assume that on the boundary of the cylinders an imedpance boundary 
condition holds, see ([5]) below. Let a — 0.5diamZ) m . Our basic assumptions are 

fca< 1, (1) 

where k is the wave number in the region exterior to the union of the cylinders, 
and 

jV(A) = — !— / N(x)dx[l + oil)], a-»Q, (2) 
lira J A 

where 7V(A) = y\ eA 1 is the number of the cylinders in an arbitrary open 
subset of the plane P, N(x) > is a continuous function, which can be chosen 
as we wish, and lixa is the arclength of a circle of radius a. The points x m 
are distributed in an arbitrary large but fixed bounded domain on the plane P. 
We denote by f2 the union of domains D m , by £1' its complement in P, and by 
D' the complement of D in P. The complement in R 3 of the union C of the 
cylinders C m we denote by C . 

The EM wave scattering problem consists of finding the solution to Maxwell's 
equations 



in C", such that 



V x H = -iueE, 
E t = ([n,H] on dC, 



(3) 
(4) 

(5) 



where dC is the union of the surfaces of the cylinders C m , E t is the tangential 
component of E on the boundary of C, n is the unit normal to dC directed out 
of the cylinders, /i and e are constants in C, lu is the frequency, k 2 = uj 2 eii, k 
is the wave number. Denote by n 2 — e/i, so k 2 — uj 2 n 2 . The solution to (J3j)- 
must have the following form 

E(x) = E (x) +v{x), x = (x-l,X2,x 3 ) = (x,y,z) = (x,z), (6) 
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where Eq(x) is the incident field, and v is the scattered field satisfying the 
radiation condition 

VF^-ifc«)=o(l), r = {x\+xlfl\ (7) 

and we assume that 

Eq(x) = k- x e iKy+ik3Z {-k^e 2 + /se 3 ), « 2 + k 2 = k 2 , (8) 

{ e j}i 3 — 1) 2, 3, are the unit vectors along the Cartesian coordinate axes x, y, z. 
We consider EM waves with H3 := = 0, i.e., E-waves, or TH waves, 

v— V x E 

E=} Eie-j, H = H x ex + H 2 e 2 = — . 9 

^— ' " noil 

i=i 

One can prove (see Appendix) that the components of E can be expressed by 
the formulas: 



%U Xj e*>*, j = 1,2, E 3 = Ue ik **, U=^< 
k z k 



where u Xj := J^-, u = u(x, y) solves the problem 

(A 2 + n 2 )u = in fl' (11) 

(u n + i£u)\ dn = 0, % := Vu-n, £ := (12) 

■w, (13) 



«a 



and w satisfies the radiation condition ([7]). 

The unique solution to ©-© is given by the formulas: 

£i = ^C4e lfe3Z , E 2 = l -^U y e ik3Z , E 3 = Ue^ z , (14) 

H 1 = -^U y e lk3Z , H 2 = — k ^—U x e ik * z , H 3 = 0, (15) 

where U x := U y is defined similarly, and u — u(x) — u(x,y) solves scalar 
two-dimensional problem (jll[) -(j!3 p . These formulas are derived in the Appendix 
for convenience of the reader. 

Problem (TTTjl - fll?)) has a unique solution (see, e.g., [I]) provided that Re( > 0, 
or, equivalently, that Im^ > 0. This corresponds to the assumption that the 
material inside the cylinders is passive, that is, the energy absorption is non- 
negative. Our goal is to derive an asymptotic formula for this solution as a — > 0. 
Our results include formulas for the solution to the scattering problem, deriva- 
tion of the equation for the effective field in the medium obtained by embedding 



3 



many thin perfectly conducting cylinders, and a formula for the refraction co- 
efficient in this limiting medium. This formula shows that by choosing suitable 
distribution of the cylinders, one can change the refraction coefficient, one can 
make it smaller than the original one. 
The paper is organized as follows. 

In Section 2 we derive an asymptotic formula for the solution to (|TTj) - (fl"3"f 
when M = 1, i.e., for scattering by one cylinder. 

In Section 3 we derive a linear algebraic system for finding some numbers 
that define the solution to problem (TTTT) - (p~5|) with M > 1. Also in Section 3 we 
derive an integral equation for the effective (self-consistent) field in the medium 
with M(a) — > oo cylinders as a — > 0. At the end of Section 3 these results are 
applied to the problem of changing the refraction coefficient of a given material 
by embedding many thin perfectly conducting cylinders into it. 

In Section 4 conclusions are formulated. 

In Appendix formulas (fT4")) - (fl"5')) are derived. 

2 EM wave scattering by one thin perfectly con- 
ducting cylinder 

Consider problem (fTTj) -(fT"3 |) with f2 = D\, VI 1 being the complement to D\ in 
R 2 . Assume for simplicity that D\ is a circle x\ + x\ < a 2 . 
Let us look for a solution of the form 

u=e iK y+ f g(x,t)a(t)dt, g(x,t) := jH^( K \x - t\), (16) 
JSi 4 

where Si is the boundary of D±, Hq is the Hankel function of order 1, with 
index 0, and a is to be found from the boundary condition (jl2[) . It is known 
that as r — > one has 

g{Kr) = a(n) + ^- In - + o(l), a := a(n) := ~ + ^- In -, (17) 

and 

9(kr) = |yj^-t)(l + (l)), r^oo. (18) 

Thus, 

u = u + g(x, 0)Q + o(—=), r -} oo; Q := a(t)dt. (19) 

The condition r — > oo is satisfied numerically if r >> a. Consequently, it is 
sufficient to find one number Q in order so solve the scattering problem fjl 1 [) - 
(fl3|) for one thin cylinder. The function (fl6|) satisfies equations (fTTj) and ([T3|) for 
any tr, and if a is such that function (|16[) satisfies boundary condition (|12p . then 
u solves problem (ITT |) - (fT3|) . We assume a sufficiently smooth (Holder-continuous 
is sufficient). 
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The solution to problem (|TT]) - (fT5)) is known to be unique (see, e.g., [1]). 
Boundary condition (|12l) yields 

- u 0n (s) - i£u = i£aQ + i£ / g (s, t)a(t)dt + (Act - a)/2, (20) 

where Aa := J Si d9 g„ dt, the formula for the limiting value on Si of the 
exterior normal derivative of the simple layer potential J s go(x,t)a(t)dt was 
used, and 

«o(s) :=e lKS >, seSi; g (s, t) = ^- In — , r st :=\s-t\. (21) 

27r r s4 

If ka < 1, fc 2 = k 2 + fc 2 , then 

lto(s) = 1 + 0(/«a), U( )n — %KU2 + O(na). (22) 

Equation (f2"U|) is uniquely solvable for a if a is sufficiently small (see [S]). 
We are interested in finding asymptotic formula for Q as a — > 0, because 
u(x) in (|16p can be well approximated in the region \x\ 3> a by the formula 

u(i) = u (x) +g(x,0)Q + o(l), a -> 0. (23) 

To find asymptotic of Q as a — » 0, let us integrate equation (|2"0)) over 5i, keep 
the main terms of the asymptotic as a — > 0, take into account that 

/ dtN 2 (t) = 0, / g Q (s,t)ds = 0(a\ log a\) a -> 0, 
use formulas (|2"2"j) . and obtain 

Q = itu (xi)\Si\, (24) 

where xi is a point inside Di, \Si\ is the length of Si, \Si\ = 2na if Si is the 
circle \x\ = a, and r st = \s — t\. The reader can find the proof of the estimate 
f s go(s,t)ds = 0(a\ loga|) as a — > 0, where s,t € Si, in [14]. From formulas 
(pM)) and ([in]) the asymptotic solution to the scattering problem (|TT|l - (fT3^) in the 
case of one circular cylinder of radius a, as a — > 0, is 

u{x) ~ Mo(i) + i2ita^g{x, O)ito(xi), a — > 0, |i| > a. (25) 

Electromagnetic wave, scattered by the single cylinder, is calculated by formulas 
(fT^|) - ([T5|) in which u = u(x) := u(xi,x 2 ) is given by formula (|2"5j) . 

3 Wave scattering by many thin cylinders 

Problem (fTTj) - ([T3l) should be solved when fi is a union of many small domains 
D m , f2 = U^f =1 £> m . We assume that D m is a circle of radius a centered at the 

point X m - 
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Let us look for u of the form 

M . 

u[x) = ug(x) + ^ / g(x,t)a m (t)dt. (26) 

We assume that the points distributed in a bounded domain D on the 

plane P = xoy by formula ©. The held iio(i) is the same as in Section 2, 
uq(x) = e lfiy , and Green's function g is the same as in formulas (IT6|) - (IT8l) . It 
follows from © that M = M(a) = 0(ln ±). We dehne the effective held, acting 
on the Dj by the formula 

?i e = v,V' = u(x) — I g(x,t)aj(t)dt, \x — Xj \ > a, (27) 
which can also be written as 

M r 

u e (x) = u (x) + y2 / g(x,t)a m (t)dt. 

, ' / ■ J S m 

We assume that the distance d — d(a) between neighboring cylinders is much 
greater than a: 

d > a, lim = 0. (28) 
a^o d(a) 

Let us rewrite (f2"6")) as 

M M 



U = Uq + 



g(x, x m )Q m + / - 3(i,i m )]cr m (i)dt, (29) 

, 1 ™ 1 « Sm 



m— 1 m— 1 

where 

(t)dt. (30) 

As a — > 0, the second sum in (|2T))) (let us denote it £ 2 ) is negligible compared 
with the hrst sum in (|2"9"|) . denoted Si, 

|S 2 |«|S 1 |, a^O. (31) 

The idea of the proof of this is similar to the one given in [BJ for a problem in 
M 3 . let us sketch this proof. 
Let us check that 

\g(x,x m )Q m \ > I / [g(x,t) - g(x,x m )]a m (t)dt\, a -> 0. (32) 



If — i m | 3> 1, and fc > is fixed then 

\g(x,x m )\ = 0{ _ 1 /2 ), |flf(x,i) -5(5,x m )| = 0( _ f /2 ), 

|X X m | \X X m | 



6 



and Q m ^ 0, so estimate (f3"2j) holds. 
If 

\x — x m \ ~ d > a, 

then 

|p(x,x m )| = 0(- T ), |g(M) -5(5,£C m )| = O(-), 
m- a 

a 

as follows from the asymptotics of Hy(r) = O(lni) as r -> 0, and from the 

formulas = -#i(r) = 0(±) as r ->■ 0. Thus, (J35) holds for |x - x m | > 

d > a. 

Consequently, the scattering problem is reduced to finding the numbers Q m , 
1 < m < M. 

Let us estimate Q m asymptotically, as a — > 0. To do this, we use the exact 
boundary condition on S m and an arguments similar to the one given in the case 
of wave scattering by one cylinder. The role of the incident field uq is played 
now by the effective field u e . The result is a formula, similar to (fM|) : 

Qj — i2ira^u e (xj), a — >• 0. (33) 

Formula, similar to (|25|) . is 

M 

u{x) ~ Uq{x) + i2ira£ g(x, x m )u e (x m ) , a — > 0. (34) 

m— 1 

The numbers u e (x m ), 1 < m < M, in (I34[) are not known. Setting x = Xj in 
(|34|). neglecting o(l) term, and using the definition (|27|) of the effective field, 
one gets a linear algebraic system for finding numbers u e (x m ): 

u e (xj) = ugjxj) + i2ira4 }^ g(xj,x m )u e (x m ), 1 < j < M. (35) 

This system can be solved numerically if the number M is not very large, say 
M < O(10 3 ). 

If M is very large, M = M{a) — > oo, a — >■ 0, then we derive a linear integral 
equation for the limiting effective field in the medium obtained by embedding 
many cylinders. 

Passing to the limit a — > in system (|3"51 is done as in [T3]. Consider a 
partition of the domain D into a union of P small squares A p , of size b = b(a), 
b ^> d ^> a. For example, one may choose b — 0(a 1 / 4 ), d = 0(a 1 ^ 2 ), so that 
there are many discs D m in the square A p . We assume that squares A p and 
A q do not have common interior points if p ^ q. Let y p be the center of A p . 
One can also choose as y p any point x m in a domain D m C A p . Since u e is 
a continuous function, one may approximate u e (x m ) by u e (y p ), provided that 
x m C A p . The error of this approximation is o(l) as a — ¥ 0. Let us rewrite the 
sum in (|35|) as follows: 

p 

27ra g(xj,x m )u e (x m ) = ^ 9{X],y P ) u e{y P )^a 2J 1= (36) 

m^j P=l CCrnGAp 
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and use formula ([5]) in the form 

27ra 1 = N (yp)W\[ 1 + oi 1 )}' a ^°- ( 37 ) 



i m eA„ 



Here |A p | is the volume of the square A p . 
From (j3"oT) and (1571) one obtains: 



p 



2™ ^2 g(xj,x m )u e (x m ) = ^2 9{x ,y p )N{y p )u e {y p )W\[l +o(l)]. (38) 

m^j _ p=i 

The sum in the right-hand side of formula (|38|) is the Riemannian sum for the 
integral 



P =i 



bm^g(xj,yp)N(yp)ue(y p )\A p \= / g{x,y)N{y)u{y)dy, u(x) = Kmu e (x). 



(39) 

Therefore, system (1351) in the limit a — >■ yields the integral equation for the 
limiting effective field 

u(x) = u (x) + i£ / g(x,y)N(y)u(y)dy. (40) 



One obtains system (|35|) if one solves equation (|40|) by a collocation method. 
Convergence of this method to the unique solution of equation (j40|) is proved in 
[TO] . Existence and uniqueness of the solution to equation (|40]) are proved as in 
[5], where a three-dimensional analog of this equation was studied. 

One has (A 2 + n 2 )g(x, y) — —S(x — y). Using this relation and applying the 
operator A 2 + k 2 to equation (|40|) yields the following differential equation for 
u(x): 

A 2 u(x) + k, 2 u(x) + i^N(x)u(x) =0 x e K 2 . (41) 

This is a Schrodinger-type equation, and u(x) is its scattering solution corre- 
sponding to the incident wave Uq = e lKy . 

Let us assume that N(x) = N is a constant. One concludes from (l4"Tj) that 
the limiting medium, obtained by embedding many perfectly conducting circular 
cylinders, has new parameter k 2 n := k 2 + i£N. This means that k 2 = k 2 + k 2 
is replaced by k := k 2 + i£N. The quantity k$ is not changed. One has 
k 2 = uj 2 n 2 , k 2 = lo 2 tlq. Consequently, n 2 /n 2 , — [k 2 + i^N)/k 2 . Therefore, the 
new refraction coefficient n 2 is 

n 2 = n 2 {l + l iNk- 2 ), i = ^f- (42) 

Since the number N > and the impedance C & r e at our disposal, equation 
(l4l?l) shows that choosing suitable N one can create a medium with a desired 
refraction coefficient. 

In practice one does not go to the limit a — > 0, but chooses a sufficiently 
small a. As a result, one obtains a medium with a refraction coefficient n 2 , 
which differs from (|42l) a little, lim^o n 2 = n 2 . 
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4 Conclusions 



Asymptotic, as a — > 0, solution is given of the EM wave scattering problem by 
many perfectly conducting parallel cylinders of radius a. The equation for the 
effective field in the limiting medium obtained when a — > and the distribution 
of the embedded cylinders is given by formula The presented theory gives 
formula (|42p for the refraction coefficient in the limiting medium. This formula 
shows how the distribution of the cylinders influences the refraction coefficient. 

5 Appendix 

Let us derive formulas (IT4l - (IT5"1) . Look for the solution to ©-CP of the form: 

E 1 = e ik ° z E 1 (x,y), E 2 = e lk * z E 2 (x,y), E 3 = e^ z u{x, y), (43) 

H 1 = e ik3Z H 1 (x,y), H 2 = e iksZ H 2 {x, y), H 3 =0, (44) 
where k 3 = const. Equation j3|) yields 

u y - ik 3 E 2 = iujfiHi, —u x + ik 3 Ei = ilouH 2 , E 2 , x = E\ tV , (45) 
where, e.g., Ej, x := Equation yields 

ik 3 H 2 — iueE\, ik 3 H\ — —iujeE 2 , H 2 ^ x — H\^ y — —iuieu. (46) 
Excluding Hj, j — 1,2, from (|4"5)) and using (|4*S)) . one gets 

Ei = ^-u x , E 2 = l -^-u y , E 3 = u, (47) 

~ Ull ~ k U f_ ~ . . 

Hi = - \, H 2 = -- x -u Xl H 3 = 0. (48) 

Since Ej = E e lkiZ and Hj = Hje ik3Z , formulas ([T3 |) -(fI5 ]l follow immediately 
from pTjl-pg]). 



9 



References 

[1] L. Landau, E. Lifshitz, Electrodynamics of continuous media, Pegamon 
Press, London, 1984. 

[2] P. Martin, Multiple scattering, Cambridge Univ. Press, Cambridge, 2006. 

[3] C.Mei, B. Vernescu, Homogenization methods for multiscale mechanics, 
World Sci., New Jersey, 2010. 

[4] A.G.Ramm, Scattering by obstacles, D.Reidel, Dordrecht, 1986. 

[5] A.G.Ramm, Wave scattering by small bodies of arbitrary shapes, World 
Sci. Publishers, Singapore, 2005. 

[6] A.G.Ramm, Many-body wave scattering by small bodies and applications, 
J. Math. Phys., 48, N10, (2007), 103511. 

[7] A.G.Ramm, Distribution of particles which produces a "smart" material, 
Jour. Stat. Phys., 127, N5, (2007), 915-934. 

[8] A.G.Ramm, Electromagnetic wave scattering by small bodies, Phys. Lett. 
A, 372/23, (2008), 4298-4306. 

[9] A.G.Ramm, Wave scattering by many small particles embedded in a 
medium, Phys. Lett. A, 372/17, (2008), 3064-3070. 

[10] A.G.Ramm, A collocation method for solving integral equations, Internat. 
Journ. Comp. Sci and Math., 3, N2, (2009), 222-228. 

[11] A.G.Ramm, A method for creating materials with a desired refraction co- 
efficient, Internat. Journ. Mod. Phys B, 24, 27, (2010), 5261-5268. 

[12] A.G.Ramm, Materials with a desired refraction coefficient can be created 
by embedding small particles into the given material, International Journal 
of Structural Changes in Solids (IJSCS), 2, N2, (2010), 17-23. 

[13] A.G.Ramm, Wave scattering by many small bodies and creating materials 
with a desired refraction coefficient, Afrika Matematika, 22, Nl, (2011), 
33-55. 

[14] A.G.Ramm, Scattering of electromagnetic waves by many this cylinders, 
Results in Physics, 1, (2011), 13-16. 



10 



